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SUMMARY 

The purpose of the present investigation was two-fold. The first 
objective was to obtain a simplified calculation technique that would 
replace the complicated automatic computer calculations. The second 
objective was to obtain experimental data for single droplets vaporizing 
at pressures greater than 1 atmosphere and to check computed histories 
against these data. 

A method of calculation using reduced curves and nomograms is ex- 
plained. Graphs and nomograms that can be used to calculate histories 
of n-decane droplets are presented. Similar graphs may be constructed 
for any pure fuel with known physical properties, as described in this 
report. The accuracy of the calculation technique 1 b shown to be within 
the agreement between experiment and theory. The method of application 
of such simplified calculations to fuel sprays Is discussed, and the 
limitations of the calculations are pointed out. 

Steady-state temperatures of hexane, decane, and hexadecane were 
obtained experimentally for various temperatures and pressures from 1 to 
5 atmospheres. The experimental data agreed with the calculated values 
to within ±10° R. A few experimental unsteady- state histories for n- 
decane at 2 and 4 atmospheres pressure are presented. These histories 
show no great deviation from the calculated histories. 

INTRODUCTION 

An understanding of the processes that occur in the atomization, 
acceleration, and vaporization of fuel droplets in a spray Is important 
to the design of combustors. Atomization results in breaking the mass of 
the liquid fuel into a spray of tiny spherical droplets and thus tremen- 
dously increases the fuel surface area exposed to the air. The total 
mass vaporized from the spray may then be obtained from a study of the 
vaporization of each droplet, or a number of Individual droplets, in the 
spray. Since rapid vaporization is essential to ensure a rapid rate of 
heat release, an understanding of the phenomena occurring in the spray 
Is important. In addition. If the largest droplets do not have time to 
vaporize sufficiently, they may leave the combustion zone in liquid 
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form; this represents a loss in efficiency charged against the combustor. 
Thus, it is of interest to he able to estimate rapidly the histories of 
single droplets in order to approximate spray histories. 

Earlier investigations 'at the University of Wisconsin were concerned 
with the problem of predicting the temperature and mass histories of 
single droplets. A theoretical investigation (ref. l) showed the effects 
of the unsteady- state portion of the vaporization histories to be of con- 
siderable importance. The theoretical study was followed by an experi- 
mental study (refs. 2 and 3} primarily concerned with vaporization under 
atmospheric pressure. Liquid-fuel droplets were hung on a spherically 
shaped thermocouple junction. Heated air was then suddenly caused to 
flow over the drop. The temperature and size of the drop were recorded 
as a function of time. The experimental results were compared with com- 
puted histories; theoretical and experimental mass curves agreed to 
within 20 percent. 

The investigation reported herein is concerned with two fundamental 
problems. The first problem is to compare the theoretical predictions 
with experimental data taken at pressures higher than 1 atmosphere. The 
second problem is to evolve a simplified calculation technique that can 
be used to obtain droplet histories easily. The equations for heat 
transfer, mass transfer, and aerodynamic drag may thus be solved simul- 
taneously to result in a reasonably simple and accurate method of pre- 
dicting spray vaporization behavior. It should be noted that such cal- 
culations give no information concerning the possible chemical interaction 
between fuel and air, nor do they take into account the effect on the 
vaporization of any fuel vapor that may be mixed with the air (refs. 3 
and 4). These latter problems are. however, beyond the scope of the 
study of single -droplet vaporization. 

This investigation was carried out at the University of Wisconsin 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. Special acknowledgment is given to 
Mr. Wolfgang Rockenhauser for his help in obtaining some of the pressure 
data presented herein. 
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SYMBOLS 

surface area of liquid drop 
lumped physical parameter (table I) 
specific heat of air at constant pressure 


Cp,a 

c p,f 


specific heat of fuel vapor at constant pressure 
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°P,L 

°P,n 

D 

F(X) 

S(T W ) 

h 

K(T W ) 

K a 

K f 

Km 

M 

M a 

Mf 

m 

Nu' 

p f,L 

Fr 

P T 

Qv 


*i 

*2 


specific heat of liquid fuel 

specific heat of air -vapor mixture 

lumped physical parameter (table I) 

diffusion coefficient of air -vapor system 

integral function of X 

lumped physical parameter 

coefficient of heat transfer 

lumped physical parameter 

thermal conductivity of air 

■thermal conductivity of fuel vapor 

average -thermal conductivity in air-vapor mixture 

mean molecular weight 

molecular weight of air 

molecular weight of fuel 

mass of drop at time 9 

Nusselt number for mass transfer 

liquid vapor pressure 

Frandtl number 

total pressure 

heat received at drop surface 

heat-transfer rate 

heat of Vaporization 

first dummy scale on nomogram 

second dummy scale on 9 £ nomogram 
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Re 

Rf 


Sc 


Tl 

^m 

V 

X 

z 

e 

e t 

0 t 

p a 

p f 

p m 

Pa 

p L 

Pm 

X 

CD 


Reynolds number 

gas constant for fuel 

radius 

Schmidt number 
dummy variable 
air temperature 
liquid temperature 
mean temperature 
steady-state liquid temperature 
air velocity relative to drop 
(Nu> - 2) 

physical parameter (appendix D) 
time, sec 

total droplet vaporization time 

total droplet vaporization time for drop initially at 
steady-state conditions 

viscosity of air 

viscosity of fuel vapor 

mean viscosity of air-vapor mixture 

density of air 

density of liquid ... 

mean density 

latent heat of vaporization 
mass vaporization rate from drop 
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Subscripts: 
a air 

d decane 

f fuel 

L liquid 

m mean value 

u unsteady- state condition 

w steady- state condition 

0,1,2, . . . 0, initial value; 1,2, . . . following values or increments 

Superscript: 

* reduced value of variable, a dimensionless quantity 


EXPERIMENTAL INVESTIGATION 
Experimental Apparatus 

The basic components of the experimental apparatus described in 
reference 3 were retained for the experiments at elevated pressure (fig. 
l). A description of the optical system, air heater, and temperature 
recorder is found in reference 3 and, for the most part, is not repeated 
herein. 

The major addition to the previously used apparatus was a pressure 
bomb or flow chamber, which was fastened to the air heater and enclosed 
the air nozzle located at the bottom of this chamber. A throttled port 
at the top of the chamber allowed the outward airflow to be regulated. 

By proper adjustment of the line pressure. Inlet flow, and exit flow, a 
range of air pressures and flew velocities could be achieved. The ther- 
mocouple (fig. 2) extended down into the chamber from the top and could 
be adjusted so that the thermocouple bead was located immediately above 
and at the center of the air nozzle. The airflow over the thermocouple 
bead was controlled by an externally triggered deflecting plate. 

The static air pressure in the chamber was measured with three cali- 
brated pressure gages whose maximum, ranges were 15, 30, and 60 pounds 
per square inch gage, respectively. 
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The air velocity at the liquid droplet was measured with an Illinois 
Testing Laboratories Velometer having a range of 0 to 600 feet per min- 
ute. Since this velometer was designed to operate at atmospheric pres- 
sure, the velcsneter case and tubes had to be maintained at the same 
pressure as that of the air stream whose velocity was being measured. 

Since the velometer could not withstand prolonged periods in a hot atmos- 
phere, it was necessary to enclose the velometer in a separate housing. 

The static pressure within this housing was made equal to that in the 
chamber by means of an interconnecting tube. A pressure gage mounted on 
the velometer housing was used to compare housing and chamber pressures. 
The metal probes of the velometer were introduced into the chamber through 
a bushing and could be slid over the nozzle or made flush with the inside 
wall of the chamber. The interconnecting tube, probes, and brackets were 
all water-cooled so that the velometer housing remained cool despite the 
elevated temperature of the chamber. 

The optical system was not revised for taking pressure data except 
that the second condensing lens and the f/l6 projection lens were placed 
inside the bomb. Two 1-inch-thick glass windows made up two walls of the 
bomb and thus allowed the light beam to pass through the bomb. A small 
mirror was mounted outside the bomb in such a way as to allow the inside 
of the bomb to be observed directly. 

At 1 atmosphere the fuel droplet could be placed directly on the 
thermocouple, as described in reference 3. For data taken at elevated 
pressures a more elaborate system was required. This system is Bhown in 
figure 3. The fuel was contained in the piston device shown in the fig- 
ure. A fine and coarse thread allowed adequate control of the feed. 

Since the piston device was external to the chamber, fuel changes could 
be made without disturbing the air pressure in the chamber. It was found 
necessary to cool the fuel by means of a water jacket so that no vapori- 
zation or heating of the fuel would occur in the fuel line. Without this 
cooling, heat conducted from the bcrnb would vaporize the fuel before the 
drop could be suspended on the thermocouple bead. The hypodermic needle 
could he brought up to the thermocouple by means of the threaded disk 
shown in figure 3. The needle could be drawn back most rapidly by manu- 
ally sliding the entire fuel line and water jacket hack an inch or two. 
When changing fuels, it vas necessary to wash the entire fuel system 
carefully with acetone before introducing the new fuel. 

Numerous attempts were made to add cool air to the chamber either 
prior to or during the experiment. It was thought that such subsidiary 
air would help to cool the thermocouple and fuel prior to taking the 
history. None of these attempts was fully successful. Any air intro- 
duced at the sides or bottom of the chamber caused fluctuations in tem- 
perature and velocity at the test section. A method was devised to turn 
on the subsidiary air anfl simultaneously open an additional air release 
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port. This method prevented disturbance of the air pressure and tempera- 
ture but caused a disturbance of the velocity that lasted approximately 
30 seconds after the cool air was turned off. Since the cooling effects 
lasted only slightly longer than the velocity disturbance, this method 
was impractical. 


Procedure and Accuracy 

The first step in taking experimental data was to calibrate the 
optical system as described in reference 3. After this calibration, the 
thermocouple was placed in the bomb, and the system was brought up to the 
desired temperature and pressure. The system usually required 24 hours 
of heating to establish a constant air temperature. Before data were t ake n, 
the air velocity and pressure were measured and recorded. The pressure 
measurements had a relative error of less than 2 percent. The air tem- 
perature was measured by using the thermocouple itself. The thermocouples 
were calibrated in a hot oil bath prior to being used; this calibration 
was accurate to within £2° R. The procedure in taking data consisted of 
the following actions: Slowly bringing the hypodermic needle up to the 

bead, placing the drop on the thermocouple, pulling back the needle 
quickly, starting camera and recorders, and triggering the air blast. 
Accuracy and calibration data for the recorder, film, and temperature 
readings may be found in reference 3. 

Several runs usually were made with the same fuel to check repro- 
ducibility. The unsteady-state portion of the histories proved to be 
less reproducible than would be desired. The thermocouple became heated 
because of the continued flow of hot air in the confined space of the bomb 
prior to taking the history. This heating, together with entrapped air 
bubbles in the drop, caused the initial portion of the temperature curve 
to vary considerably from run to run. Repeated trials produced curves 
believed to be reasonably accurate in the unsteady- state portion of the 
history, but the difficulties of the experiment prevented the taking of 
data that could be considered as reliable as that presented in reference 
3. Reproducibility of steady- state temperatures was found to be good, 
with variations as small as ±3° R. 


Experimental Results 

A limited amount of acceptable data was obtained for pressures 
greater than 1 atmosphere. Several trends, however, were experimentally 
established. As the pressure is increased, the unsteady- state becomes 
even more important; that is, the ratio of the time to attain the steady- 
state temperature to the time to completely vaporize the drop increases. 
Theoretically, this ratio approaches a limit of 1 as the pressure is 
increased. Secondly, the steady- state temperatures increase with in- 
creasing pressure. This increase tends to become less at higher pres- 
sures, so that the difference between unsteady-state temperatures at 1 
and 2 atmospheres is much greater than the difference at 4 anti 5 
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atmospheres. Thirdly, the higher steady- state temperatures may in some 
cases cause the total vaporization time to decrease with increasing 
pressure. 

Figures 4 and 5 show experimental temperature and mass histories at 
2 and 4 atmospheres of pressure. Figure 6 gives experimental steady- 
state temperatures for various pressures against air temperature. At 
higher air temperatures, heat transfer from the thermocouple wires caused 
the liquid temperatures to he erroneously high. This heat-transfer prob- 
lem was increased at higher pressures because the thermocouple itself was 
heated more in the bomb than when in open air. 


Comparisons of Experimental and Calculated Results 

Figure 6 shows the agreement between calculated and experimental 
steady-state temperatures. This agreement is within ±10° R as compared 
with +4° for 1 atmosphere of pressure. About ±5° of this difference may 
be accounted for by experimental inaccuracies. Figures 4 and 5 show 
both computed and experimental histories at 2 and 4 atmospheres. The 
same error trends that were discussed in reference 3 are presentj that 
is, the experimental and mass-transfer histories were within 20 percent 
for all cases studied. With the exception of steady-state temperature 
agreement, the agreement between experimental and calculated histories 
was similar to that reported in reference 3. No other significant error 
trend was found as the pressure was increased to 5 atmospheres 
absolute. ' 


SIMPLIFIED CALCULATION TECHNIQUE 

Original computations of droplet vaporization histories, as found 
in references 1 to 3, were performed with the aid of either desk Calcu- 
lators or automatic high-speed computers. The desk calculations were 
prohibitively lengthy. The time required to compute histories with 
sequenced computers varied considerably with the type of machine avail- 
able. The IBM 650 with subroutined floating decimal can calculate a 
history in approximately 10 minuteB. Other machines may require as much 
as l/2 hour or as little as 2 minutes to calculate a history. 

While the high-speed computer is an excellent tool for calculation, 
it has a number of disadvantages: Computer programing may be time con- 

suming, computer time may be difficult to obtain, the cost of operation 
is high, and lastly a computer of adequate size simply may not be avail- 
able. Therefore, an attempt was made to replace the computer calcula- 
tions with a simplified calculation technique. A number of techniques 
were tried with little successj these attempts are reported in appendix 
A. The method of reduced- variable graphs, which proved most successful, 
is described herein. 
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Construction of Plots of Reduced Variables 

The mass-transfer histories given in reference 3 are quite similar 
in appearance, the characteristics curve being "S” shaped. Deviation 
from this shape occurs only when the initial slope is not approximately 
zero. Such deviations may "be eliminated by restricting the initial liquid 
temperature to same appropriately low value. The exact value of the liq- 
uid temperature that will .satisfy this condition depends on the fuel and 
air temperature range. For example, the ma ximum initial liquid tempera- 
ture for n-decane with air temperatures of 800° to 1600° R that will give 
an initial slope of approximately zero is 550° R. This maximum, 'tempera- 
ture is c all ed the zero- slope temperature. 

With the initial temperature fixed at the zero-slope value, all 
histories for a particular fuel are similar in appearance except for the 
time scale. For a particular fuel under the conditions studied, division 
of the time scale by the total vaporization time gives a family of reduced- 
mass curves that all lie within 20 percent of each other. Xf the air 
temperature is held constant and the radius and/or velocity varied, the 
reduced curves fall very close to each other. Variation of the air tem- 
perature thus causes most of the 20 "percent variation mentioned. 

Examination of the mass curves shows that the inflection point cor- 
responds to the end of the unsteady state. The proportion of time spent 
in the unsteady state is known to be influenced by the air temperature. 
Fortunately, the reduced-mass curve is nearly a straight line in the 
neighborhood of the inflection point and thus minimizes the effects of 
air temperature. A single curve thus can be made to represent all such 
reduced curves to within the accuracy warranted by the agreement 
between IBM calculations and experiment. Such a reduced curve for n- 
decane is shown in figure 7(a). This curve was constructed from computed 
mass-transfer curves. When the reduced curve was drawn, points farther 
to the right were favored so as to fit the experimental trends better, 
as discussed in reference 3. 

Computations for initial liquid temperatures lower than the zero- 
slope temperature offer no difficulty. To a good approximation the equa- 
tion for heat transfer with no mass transfer applies to these low liquid 
temperatures. Appendix A gives the equation and solution for this case. 

Computations for initial liquid temperatures higher than the zero- 
slope temperature require the use of only part of the reduced-mas s- 
transfer curve. The portion of the reduced-mas s-transfer curve to the 
right of any value of reduced time corresponds to a mass history with 
initial temperature greater than the zero-slope temperature. If the liq- 
uid temperature at each reduced time is known, the reduced-mas s-transfer 
history may be obtained by merely taking that part of the history to the 
right of the point corresponding to the desired initial temperature. New 
reduced axes that run from 0 to 1 must then be constructed so that the 
history will begin at zero time and zero mass transferred. 
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A plot of liquid temperature against time gives a single curve for 
any one fuel, pressure, and air temperature; that is, the curve is inde- 
pendent of size and velocity. Families of curves must "be drawn for var- 
ious pressures and temperatures. Figure 7(b) gives such curves for n- 
decane at 1 atmosphere of pressure. Again, t^ e curves were constructed 
from computed data. A reduced-temperature scale does not bring these 
curves for various air temperatures together. 

The total vaporization time may be obtained from steady-state vapor- 
ization time and the reduced-mass curve, as explained in the following 
section. As a computational aid, the graphs and nomogram presented in 
figures 8 to 10 were constructed from equation (Al) and the physical 
properties given in appendix IV of reference 3. Appendix B gives 
details concerning the construction of figure 10. 

Graphs and nomograms similar to those given here may be constructed 
for any pure hydrocarbon fuel. Hexane and hexadecane, for example, were 
reduced in the same manner as decane. Unfortunately, the shape of the 
reduced-mass curve varies with the fuel. Ho single curve that would 
represent all fuels could be found. When the curves for yarious fuels 
are plotted so as to agree during the initial portion, they diverge at 
higher values of the reduced time 6*. 

If it is desired to calculate a mass history for a fuel "x" similar 
to decane (i. e. , the steady-state temperatures are within approximately 
50° R of each other for several air temperatures), the reduced-mass 
curve for decane maybe corrected as follows: Use the reduced curves for 

decane with the initial liquid temperature computed from the equation 

T L,0,d = ®L,0,x + ^ T w,d " T w,x^ 

where T-^ q ^ is the initial liquid temperature to be used for decane, 

Tl © x ^■ he initial liquid temperature of fuel M x, " T w ^ is the 
steady-state temperature for decane, and T w x is the steady-state tem- 
perature for fuel "x. " The construction of the reduced curve follows the 
procedure given in the following section. The total vaporization time 
should be calculated from data for fuel "x. " It should be remembered 
that this correction gives only an approximation to the correct reduced 
curve and should be used only for fuels similar to decane. 

The reduced-mass curve for decane at 1 atmosphere pressure may be 
used with but little error for higher pressures up to at least several 
atmospheres. The liquid-temperature curves will not be the same for 
different pressures. 
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Calculation of Mass and Temperature Histories Using ELots of 
Reduced Variables for Droplets Vaporizing with. 


Constant Relative Air Velocity 


Figure 7(a) is a plot that represents all mass- transfer curves 
(within the accuracy given on p. 14) for the following range of 
conditions : 


Initial liquid temperature, °R 

Air temperature, °R 

Air velocity, in. /sec .... 

Initial radius, in 

Fuel 

Total pressure, atm abs . . . 


550 

. 800 < Tb s 1600 

. . 10SV < 1000 
4X10“ 4 < r Q £ 10“ 2 

n-decane 

. . . . 1 < P T < 4 


Figure 7(b) is a plot of liquid temperature against reduced time for 
various air temperatures. Each curve is for 1 atmosphere of total pres- 
sure and the given range in radius and velocity. 


Two steps are involved in the calculation of mass histories from the 
charts. Step 1 is the construction of a reduced-mass curve appropriate 
to the initial liquid temperature of interest. Step 2 is the conversion 
of the reduced curve to the desired mass (against time) history. The 
procedures to be. used in different canes are indicated as follows: 


Step 1; Construction of reduced-mass curves 

(a) T Lj0 = 550° R 

Drops with initial liquid temperature of 550° R have the reduced- 
mass-transfer curve given by figure 7(a). 

(b) Tj^o > 550° R 

If the initial liquid temperature is greater than 550° R, only a 
part of the reduced -mass curve given in figure 7(a) can be used. The 
construction of the new reduced-mass curve for this case proceeds as 
follows: 


(1) Select from figure 7(b) the curve corresponding to the air tem- 
perature of interest. 

(2) Read the value of 9* given by this curve at the temperature 
Tj^qJ call this value of the reduced time ©*. 
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(3) Find the point corresponding to 0* on the curve of figure 7(a): 

call, this point (m*, 0*}. 0 

(4) Let the point (mj^, 0^) he the origin of a nev set of axesj these 

axes should he constructed parallel to the original m*, 0* axes. Con- 
struct scales on each of these nev axes that run from 0 to 1. 0. The 
value 1 should coincide with the value 1 on the old m*, 0* axes. Fig- 
ure 11 is a schematic example of such a construction. The nev axes and 
the portion of the reduced-mass curve to the right of the nev origin 
(m*, 0*) form the reduced-mass curve sought. The uniform scales (fig. 

12) can he used in constructing the nev axes scale, or a straight edge 
and the geometry of similar triangles can he used in place of the unif orm 
scales. 

(c) T-^o < 550° R 

If the initial liquid temperature is less than 550° R, the portion 
of the mass history prior to 550° R must he computed separately. Since 
no mass is transferred during this period, the calculation is for the 
heating time only. This case proceeds as follovs : 

(1) Use equation (A2c) vith values of C and D as given in table 
I to calculate the time required to reach 550° R liquid temperature. 

Rote that this is not a reduced time. 

(2) Calculate the actual mass history from 550° R on as directed in 
step 1(a) and step 2. 

(3) Add the time to reach 550° R (step l(cl)) to the result of 
step l(c2). 

Step 2: Conversion of reduced-mass, curves to actual-mass histories 


The reduced-mass curve given hy step 1 is converted to a mass- 
transfer history hy multiplying the reduced time 0 q hy the total vapor- 
ization time e t> ■which is obtained as follovs: 

Pick a value of 9 * that is veil within the steady-state region 
(say 0* = 0. 8) j call this point (m^, 8^). Find the reduced radius 
r* = r w /r 0 corresponding to this point hy use of figure 13. Figure 13 
requires the knovledge of T L Q , T v , and m*, all of which are known. 

Figure 10 nov gives the vaporization time for a drop vaporizing at steady- 
state conditions with initial radius r w and velocity V. The value of 
0^ can nov he calculated from. the equation 

e t - ®t/ 1 - 

A numerical example of this procedure Is given in appendix C. 
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Construction of Mass Histories for Droplets Vaporizing with 
Changing Relative Air Velocity 

It is assumed that drag coefficients are available from, which the 
relative air velocity may he calculated (see ref. 5). The procedure to 
he followed is a stepwise trial- and- error process: 

(1) Construct the reduced-mass-transfer graph as in step 1 of the 
preceding section. 

(2) Hold all conditions except velocity constant and calculate values 
of 0-j. for various V. The calculation is performed as in step 2 of the 
preceding section. 

(3) Pick an increment of time that is some small fraction of 

the value corresponding to the given initial velocity. 

(4) Assuming the mass to he constant, calculate the velocity at time 
A0 as predicted hy the drag equations. Calculate the average velocity 
during time A 0. 

(5) From the graph prepared in step (2) of this section, read the 
value of 0 t corresponding to the average velocity calculated in step 
(4). Cal l this time 0-^,1* 

(6) Calculate A0* = A0/© t 

(7) From the mass curve of step (l) of this section, find the m* 
corresponding to A0*. This gives the first point, (m*, A0), of the 
mass history. 

(8) If the mass has changed appreciably during the time A0, the 
calculations of steps (4) to (7) may need to he repeated with a more 
accurate average mass. 

(9) Repeat steps (4) and ( 5 ) for the next time step. Step (6) now 

should read A©* = A + A 0/© t Step (7) follows as before. 

For steps other than the first, the curves for velocity and mass of 
the drop may he extrapolated as an aid in estimating the values to he 
used for the next step. 

The calculation thus proceeds step hy step to give the mass-transfer 
history. A numerical example is given in appendix C. 
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Accuracy of Method of Reduced -Variahle Charts 

Ease of calculation and speed are important factors in considering 
the relative merits of calculation techniques, but the accuracy of the 
results must not be ignored. The range of conditions covered by the 
reduced-variable plots was chosen so as to approximate engine operating 
conditions. This range exceeded the range for which experimental data 
were available. Direct comparison of mass histories obtained by use of 
the reduced-variable charts and the experimental curves was possible only 
for the largest drop size in the range considered, that is, 500 microns. 
Comparisons of experimental curves with curves obtained from IBM calcula- 
tion and from the reduced- variable charts are given in figure 14. Curves 
for calculations, assuming no unsteady state and also for the method of 
solid spheres (see appendix A) are included. In each case the use of the 
reduced-variable charts provides the best agreement with experiment. In 
order to check these charts for a wider range of conditions, comparisons 
with IBM calculations only are given in figure 15. In each case the 
agreement is good. In general the agreement is poorest at low air tem- 
peratures. Radius and velocity have but little effect on the relative 
error. While the charts agree with the IBM calculations for these condi- 
tions, it must be remembered that the theoretical calculations themselves 
are subject to the error trends given in reference 3. These trends sug- 
gest that the IBM calculations may disagree considerably with experiment 
for very small droplets. 

Figure 16 compares the results of calculations using the IBM and the 
reduced-variable charts for variable air velocity. The velocity curve 
was chosen arbitrarily and thus was known precisely. This procedure is 
permissible, since the computations were for comparison purposes only. 


Use and Limitations of Histories 

In order to use the computed histories to predict the behavior of 
an air-spray system, the following technique may be used: 

(1) A representative droplet- size distribution for the condition of 
Interest is determined or assumed. The initial liquid temperature must 
be known or assumed. The airflow pattern and relative air velocity must 
also be known or assumed. 

(2) The use of the computed histories that assume no interaction 
between droplets should be justified on the basis of the dispersion of 
the spray in question and the over -al l , fuel-air ratio. In the case of 
sprays found in jet engines, the authors believe that this assumption is 
justified and that the computations presented are adequate. If the spray 
is dense, other considerations and techniques will have to be taken into 
account (refs. 3 and 4). 
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(3) The entire spray may then be divided into a number of "size 
groups" and the percentage by -weight of each determined. Each size group 
may be represented by a single droplet size that may be chosen as some 
average of the range of its particular size group. Any error involved 
in the averaging te chn ique will be minimized if a sufficiently large num- 
ber of groups is chosen. 

(4) Each representative droplet will be treated in the manner sug- 
gested in this report. Since the droplets will be affected by aerody- 
namic drag, the procedure outlined on page 13 should be used. 

(5) The amount vaporized from each size group will then be obtained 
as a function of the distance from the nozzle or the point at which the 
spray was first formed and came into contact with the air. 

(6) The total amount of vapor formed from the spray can then be ob- 
tained by adding the properly weighted vaporization due to each size 
group. 


(7) After the initial fuel-air ratio by weight is known, a distance 
from the nozzle may be found at which an over- ail air-vapor mixture of 
combustible strength is formed and where all physical processes are 
thought to have made their contributions. A spark plug or flameholder 
may then be placed near that point to initiate or sustain combustion. 

The balance of the computed individual histories are actually meaningless 
beyond this point. 

The calculated vaporization histories do not necessarily describe 
the droplet lifetimes in their entirety. Seme or all the droplets 
present, say in a fuel-spray system, may deviate from these calculated 
histories sometime during their lifetimes in an actual combustor. The 
lifetimes of fuel droplets may he thought of as composed of two parts. 

The first is physical; this can he adequately described by the calculated 
vaporization histories obtained by using the technique presented in this 
report. The second is chemical; this includes the time required for the 
oxygen and fuel molecules to arrange themselves for combination into 
combustion products. It is quite likely that the two mentioned parts, 
physical and chemi cal, overlap to an appreciable degree. 

A second limitation has to he placed on the unrestricted -use of 
the computed histories as previously outlined. These computations are 
based on the assumption that an abundant quantity of air accompanies 
each droplet; that Is, no interaction occurs between droplets, and the 
pressure and temperature of the air are unaffected by the vaporization 
process. Thus the computed histories do not show any effect that may 
result from the close proximity of droplets. In dense parts of a 
spray the air temperature may drop appreciably. The rate of vaporization 
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also will decrease "because of the reduction of diffusional driving forces 
by the presence of fuel vapor in the air surrounding the droplets. As 
an ultimate condition, the fuel-air system will reach adiabatic equilib- 
rium, and the vaporization process ceases (refs. 3 an d 4). 

It is felt that, for the most part, fuel sprays existing in continu- 
ous flow or Jet-type combustors are so dispersed that the partial pressure 
of the fuel vapor in the air surrounding the droplets is not high enough 
to affect the diffusional driving forces appreciably. The driving force, 
therefore, will always be assumed equal to the fuel vapor pres stir e of the 
liquid fuel at the liquid temperature. The droplets are far enough re- 
moved from each other that "single droplet" computations, such as those 
described in this and previous reports, are accurate .enough to describe 
the histories of the droplets during their period of pure vaporization) 
that is, before ignition takes place. The air temperature during this 
process will not be appreciably reduced. The air pressure will remain 
uniform, since frictional losses and the pressure rise due to vaporiza- 
tion are relatively negligible. 

Caution should also be exercised in the use of "ignition delays" or 
"preflame distances" obtained by using one of these procedures. In 
actual combustors, there are a number of variables that render this pro- 
cedure comparative at best. Turbulence and backflow, for example, dis- 
turb the pattern of relative velocity between air and drops. Atomization 
uncertainties are also of concern. Heat transfer by radiation from the 
combustor walls and flame to the liquid droplets adds to the convective 
heat transfer to the droplets and is difficult to evaluate. The evalua- 
tion Just given should, therefore, be used for estimation or comparison 
purposes only. 


CONCLUSIONS 

An investigation was conducted to (l) compare the theoretical pre- 
dictions with experimental data obtained at pressures from 1 to 5 atmos- 
pheres, and (2) evolve a simplified calculation technique of reduced- 
variable graphs and nomograms so as to obtain droplet histories easily. 
The results and conclusions are stated as follows : 

1. The ratio of the unsteady- state time to the total vaporization 
time increases with increasing pressure. It is conceivable that at very 
high pressures the entire droplet lifetime may be spent in the unsteady 
state. 


2. Agreement between calculated and experimental mass-transfer 
histories was within 20 percent for all cases studied. Reasonable agree- 
ment between theory and experiment was found for pressures up to 4 
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atmospheres. Extrapolation of the trends of agreement between exper iment 
and theory suggests that the theory may be in serious disagreement with 
experiment for droplets considerably smaller than 250 microns initial 
diameter . 

3. The experimental method i6 definitely limited in scope. This 
method ceases to be practicable for very small drops or higher air ve- 
locities. Small droplets imply such short vaporization times that it iB 
impossible to plane the drops on the thermocouple before they vaporize. 

In addition, the time to develop the airstream profile may become impor- 
tant and thus raise the uncertainties of turbulence and air-velocity 
magnitude. Air velocities are limited by the effects of droplet distor- 
tion and blowoff. It is conceivable that the method may be used for 
higher air temperatures or pressures but only with considerable experi- 
mental difficulty. It may be concluded that seme other method should be 
sought if conditions greatly different from those studied herein are to 
be considered. 

4. The method of reduced-variable charts presented in this report 
is simple enough and accurate enough to be used for the calculation of 
the vaporization of sprays. The method provides as much accuracy as is 
warranted by the agreement between theory and experiment. 

5. Many uncertainties in the application of these studies to fuel 
sprays result from a lack of understanding of the detailed behavior of 
sprays. Drop-size distribution, drag coefficients, droplet interference 
effects, and the effects of chemical reaction in the boundary layer are 
some of the more important of these problems. 


University of 
Madison, 


Wisconsin, 
Wis., May 9, 


1957. 
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APPENDIX A 


OTHER METHODS OF CALCULATION 

A flow diagram of the IBM 650 program 1b given in figure 17. Values 
of T^l and r for each new time step are extrapolated from, the three 
prior points. The program automatically recycles each step until the 
assumed and calculated values agree to within a Bet difference. 

A procedure involving the same calculation technique used in the 
IBM program hut utilizing nomographic aids and graphs was devised. Be- 
cause the stepwise procedure was retained, it was estimated that at "best 
this procedure would require 0. 7 hour of errorless calculation per history. 
Such a procedure offered little advantage to the computing machine method 
and was thus abandoned. 


Appendix D presents the equations used in the calculations. The 
extreme nonlinearity and complicated coupling of these equations pro- 
hibits any hope of an analytical solution for the most general case. 
The following two Important special cases do allow solution. 


Steady-State Heat Transfer, dT^/d Q = 0 
In this case equation (Dl) can be integrated to give 

r§„F(X) 


e ; = 


f-W^- 


K 


where 




0,w 
K = 


total lifetime at steady-state condition 
initial radius at steady-state condition 


X = 


F(X) 


D^ 

inf P * - 

R f T m p l 

\ T " P f,L 

2 

1 

Nu' - 2 

0.3 ReV^c- 

Y / 

x s ds 

X i 

, 77 ?!* 


(Al) 


V 


I 
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o 

s 


in 


(See list of symbols for -undefined, quantities. ) A closed-form solution 
of the integral F is given in reference 2, which, also gives two limit- 
ing values for F. Only the limiting value for large values of X is 
appropriate to .jet engine conditions. 3h that case F approaches the 
constant value 0. 500. Graphs of g, F, and K for n-decane at 1 atmos- 
phere pressure are given in figures 8 and 9. These graphs can be used 
in conjunction with a slide rule or desk calculator to obtain 01. A 
nomographic SQlution of equation (Al) is also presented in figure 10. 

Negligible Mass Transfer, dm/d 9 = 0 
For this case equation (Dl) becomes 

S - | 5^152 ( 2 + 0. 6 ErVW/ 2 ) (A2a) 

" 2 *\%,L 

The equation is applicable to n-decane for liquid temperatures below 
550° R. For a range in liquid temperature from 460° to 550° R, the 
change in radius of a n-decane droplet due to thermal expansion is less 
than 3. 5 percent. Rearranging the equation gives 


de = 


r 2 

G + D(rV) 1 / 2 t B " T L 


(A2b) 


where 


0 - 




The properties of air are evaluated at the mean temperature, 

2^ + T t n + 550 

m _ B U 

nm 4 

A tabulation of C and D (see table I) for n-decane shows that these 
factors change only slightly over a 50° R range in T m . Thus, integrating 
gives 


- 


* C + D(rV) 1 / 2 \ B " 


T B ~ T L.Q 1 
550 i 


(A2c) 


In 



20 


NACA TN 4338 


While these two special cases are very useful, they do not directly 
give any information about the important part of the heating period. 

If the unsteady- state time Q ^ u and the steady- state temperature T^ 
are known, the following approximation to the temperature history may be 
made: 


T* s* 

ip* _ 


T L " 


T. 


w 


T- 


L,0 


5 t,u 


for 

VI 

CD 

0 t,u 

for 

CD 

IV 

®t,u 


With this approximation it is possible to obtain integrable forms of 
equation (Dl). These forms are, however, highly complicated. A further 
objection to this method is the required knowledge of 0^. u . Higher 
order approximations to the temperature histories give unintegrable forms 
and require the knowledge of an additional temperature-time point. 


Various equations were fit to the temperature and mass histories in 
an attempt to represent these curves with a universal formula. The fol- 
lowing equations give reasonable accuracy: 


«_ *w- T L -AO* 
T„ - T t n 



/ 

n \m 

1 - nr = rj- = exp 
Mq 


Vi)_ 


(A3) 


(A4) 




Values of constants A, B, n, m are needed for very good accuracy. No 
good correlation of these constants with ambient conditions could be 
found. This representation further indicated that a reduced plot of some 
kind was possible. The values of B and m indicated that a reduced 
plot should be based on a fixed initial liquid temperature. 

Another simplified calculation technique is given here for complete- 
ness. This method calculates the time to reach steady-state temperature 
assuming no mass transfer. From this point on, the steady-state vaporiza- 
tion equation is used. The calculation for the unsteady- state time thus 
ignores the effect of mass transfer on heat tr ans fer, the latent heat of 
vaporization, and the change in droplet radius. These effects may be 
large, since 40 to 50 percent of the drop may vaporize during the un- 
steady state. While the method is approximate, comparisons given in 
figure 14 show that it is a very definite improvement over the simple 
steady- state calculations. 


080S 
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APPENDIX B 


CONSTRUCTION OP 8|. NOMOGRAM 


The- e£ nomogram (fig. 10} referred to in the tody of the report 
was constructed from computed values of the functions and the physical 
properties given in reference 6. The nomogram calculates the value of 
0_k from the equations 


0 t 


r 2 

K(*w) 


P(X) 


X = (rV) 1 / 2 g ( Tw ) ( see p. 18) 

The construction with the form suggested in reference 6 is outlined 
as follows: 

(1) log (r) + log (V) = q x 

(2) 2 log (X) = q. ± + 2 log (g) 

(3) log (F) + 2 log (r) = q 2 

(4) log (0^.) = q , 2 - log (K) 


Varl- 

Range 

Scale modulus 

Plotting 

Matching point 

ahle 



modulus 


r 

4XL0" 3 to 4XL0" 3 

. 14 

14 

4X10“ 3 (hottcm base) 

V 

10 to 10 3 

14 

14 

10 (hottcm) 

*1 


14X14 _ 

14+14 



g 

3.563 to 6.000 

65 

130 

3. 563 (hottcm) 

x 


7x65 _ 6e31 94 

12.639 

7.126X10” 2 (hottcm) 



7+65 



F 

0.090 to 0.500 

40 

40 

0.090 (top base) 

r 

4x10“ 5 to 4x10 " 3 

7 

14 

4x10“ 3 (top) 

*2 


5.957 

7+40 



K 

10“ 3 to 8. 0816X10“ 5 

32 

-32 

8. 0816X10 -5 (top) 

0 t 


5. 96x32 _ c noo 

5. 022 

1.78X10" 6 (top) 


5. 96+32 




22 


NA.GA. TN 4338 


Distances "between scales: 


r X V 

arbitrary 10" 

!l r l 

10 = 5" 

28 


arbitrary 12" 

ijx 

I 12 — 1 ' 

7 + 65 


arbitrary 6. 7' 


H 6. 7 = 14" 


arbitrary 14" 

v; 

5,96 11 - 1 

5.96+32 


Since there are actually two nomograms connected through the func- 
tion F(X ) , the distance between the F and X scales is purely arbi- 
trary and does not affect the other scale distances. 

The scales from left to right as they appear on the nomogram are: 

r l> V s ' F > Y > V r 2' Q V g > K ' 

Scale construction: The F and X scales were matched from IBM 650 

computed values of F against s. The scales g, K, X, F were specially 
computed for each scale division. The scales r , V, and &{. were taken 
from logarithmic scales provided in reference 6. The scales g and 
K were marked in terms of the independent variable T w< 

Calculation proceeds as follows; 

(1) Pick the values of r^ = r 2 , V, and T w of Interest. 

(2) Draw a line from r^_ to V. 

(3) From the intersection of r-jV on q-^ draw a line to g(T w ). 

(4) Read the value of X at the intersection of q^g on X. 

(5) Find the point on the F scale that equals this value of X. 
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(6) Draw a line frcm F to r 2 * 

(7) From tiie intersection of this Fr 2 line with q. 2 draw a line 
to K(T W ). 

(8) The intersection of this q^K line with the 6^ scale gives 
the desired value 1 of 8^.. 

Construction of such a nomogram requires many hours of tedious work. 
It is therefore felt that if values of 0.J. are desired for other fuels, 
the graphs of F, g, and K should he drawn and used directly. A 
nomogram should he constructed only in cases where very extensive use is 
anticipated. 
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APEEHDIX C 


EXAMPLE OF CALCULATIONS USING VARIABLE VELOCITY AND 
CHARTS OF REDUCED VARIABLES 

As an example of the method of reduced- variable charts , the follow- 
ing initial conditions will he used: 


Air temperature , Tg, °R 1380 

Initial relative air velocity, V a ^g, in. /sec 1000 

Initial radius, rg, in 0.001 

Initial diameter, microns 50 

Initial liquid temperature, °R 600 


Ingebo ’ s drag coefficient (ref. 5) is used here: 


= 27 

„ 0.84 
Re 

££ - . 3/£a\/v£\ 

de 8 Vp l ]W^ 

Step 1 : Using only pant of the curve of figure 7(a), obtain the reduced 

curves shown in figure 18. Hie origin of the new axis is 0.065 on the 
curve of figure 7(a). Thus, only the 9* scale must he revised. For 
example, 15 percent is vaporized at 0* = 0.25 on the new scale instead 
of 0* =0.3 as on the old scale. 



Step 2 : Obtain values of 0^ for various velocities. For 9* = 0.8, 

T w = 736, m* = 0.9 calculate r * = 0. 475, ~r w = 0. 000475j 0 t = 0^/O.2 = 50^ 
Thus, ©-j-XLO^ = 4.42, 4.62 , 4.75, and 4.88 seconds; V = 1000, 900, 800, 
and 700 inches per second. Plot these values as in figure 18. 

Step 3 : Pick A0 = 5x10“ 4 . 

Step 4: Calculate the velocity at time A 6 by using the drag equations: 


V - v = — P ' & 

V I v 2 4-1.84 


y 0.84 

71. 16 Ae _ K Agjyl. 16 


where 

V 1 = 1000; Vg is the velocity at time A 0; V = (V-j_ + Vg)/2 
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p (1380) = 1. 66X10" 5 

r = 0. 001 

P l (625) = 0.02487 (estimating AS* will be * 0.1) 
v a = y a (!000) = 0. 0646 


t 

s 


For a starting value, use V P = ^ After this, extrapolation of 

1 + 

u 

the velocity curve should, he used to obtain Tfa, 

= 5X1. 66X10~ 5 (0. 0646) 84 
2. 487X10 -2 

K = 16. 6 
V = 990 

V g ■ V x - 16. 6(3000) (5XL0“ 4 ) = 1000 - 25 - 975 
V-i + Vo 

Recalculate V = 988. 

Step 5 : Qj. i = 4.47 from the graph of step 2. 

Step 6: AS* = = -^*3 = 2—. . = 0.113; 6* = 0.113 

1 0 t,l 4. 47X10" 3 1 

Step 7 : The first point of the m*, 6 curve is (0.025, 5X10 -4 ). 

Step 8 : The guesses of Tq and T^_ were very close to the calculated 
values, and thus no recalculation is necessary. Repeat, beginning with 
step 4, - to find the second point of the m*, 9 curve. 

V 2 = 975,* v 3 = ? 

p (1380) = 1. 66XL0" 5 

r = 0.0096 (examine r* graph; assume A 9* = 0. l) 
p L (680) = 0.02397 
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V a (l030) = 0.0686 

V = 963 (linear extrapolation of previously calculated points) 


7.5x1. eexfes.e) 0,84 ^" 1 
0.93X2.397 


19.6 


V 3 = 975 - 29.4 = 946 
9. a 4. 53X10“ 3 

*C y dt 

AQ* = — 5 - — Q — — = 0. 110! 8* = 0. 223 
2 -3 2 

4. 53X10 

The second point of the m*, 9 curve is (0.125, 10" 3 ). The guess of 
T^ was quite good; r was close enough. Clearly on the next step r 

will he approximately 0.00094, and so forth. Repeat to obtain third 
point, and so forth. 
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APPENDIX D 

RESUME OF EQUATIONS USED IN CALCULATIONS 


The following equations were used in the calculations of histories^ 
they are derived and discussed in references 1 and 2. They are repeated 
here for convenient reference! 



(D2) 


where 


m = 4/3(itr 3 p L ) 
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D y , \x & , p.j», K & , K^, Cp^ a , and. Cp^f a^re taken at temperature T^. 

P f p L , Cp and ~k are taken at the liquid temperature T^, •which 

is assumed uniform throughout the drop (ref. 2). The equations for the 
physical properties may he found in reference 7 . 
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TABLE I. - VALUES OF C 



Tm = 


2T B + T L,0 + 550 


(See appendix A) 


•] 


AND D PARAMETERS 

) = r2 
t C + D(rV) 1 / 2 


Tin* 

°R 

C 

1075 

1.378x10' 

1050 

1.350 

1025 

1.320 

1000 

1.290 

975 

1.260 

950 

1.230 

925 

1.200 

900 

1.170 

875 

1.140 

850 

1.110 

825 

1.080 

800 

1.050 

775 

1.020 

750 

.993 

725 

.965 ' 

700 

.934 

675 

.905 

655 

.880 


D 


1.81X10 -4 

1.81 

1.81 

1.81 

1.809 

1.80 

1.788 

1.785 

1.785 

1.785 

1.785 

1.785 

1.785 

1.785 

1.785 

1.785 

1.785 

1.785 



Thermocouple ■ 


Air supply -— □a 
100 lb/sq in. gage 


To pressure gages 
Water-cooled Jackets *\ \ 

Velometer-v \ Fzd 


Porous 

stone 

filter 


Filtrate 

collector 


f Pressure 

Pressure V J regulator 
eage>— v > f 


Critical 

flow 

orifices 



L-Ji 


.To throttling valves 
Pressure chamber 
/ /~Fuel line 







Water-cooled block 

Flow deflector and microswitch 

Solenoid trigger 

/ — To air-blast 
timing light 


Flow atraightenars 


« m> ill 

i y l 


Variac 


■ Steel 
beads 


110-Volt 

line 


Figure 1. - Experimental apparatus. 
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Exploded viev of "B" 



■Outer Insulator 
- Genlflex insulator 
•^Wire 


Figure 2. - Thermocouple construction used vlth pressure data 





.2 
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Figure 4. 
solute; 


Time , sec 

(a) n- Decane j initial diameter, 944 microns. 

- Exper ime ntal and calculated histories. Total pressure, 2 atmospheres ab ■ 
air temperature, 900° E; air Telocity, 130 inches per second. 




Percent mass transferred Liquid temperature 
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0 .5 1.0 1.5 2.0 

Time, sec 


Ob) Initial radius, 555 microns. 
Figure 5. - Concluded. 
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figure 6. - 
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Air temperature, PR 
(c) Eexadecane. 

"WH mn^ H _ — 

* *5“ ■» W ■ W W ■” ■ 
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Reduced tlo#, 

(b) Reduc«d-teap«rrtur* oar t*. 

Figure J . - Seduced maBB and temperature curves for n-decane. 
(A large vorking copy of this fig. may he obtained by using 
request card bound in back of report.) 



Steady- state temperature, R 
( a) Factor K. 


Figure 8. - Various factors plotted against steady-state 
temperature Ty for calculation of vaporization time. 
Total pressure, 1 atmosphere. 





Factor 



Factoj 
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Dr*» 21 m fro* ootlo *** to 'E*, pivot ohoro AE orotaoo "B" to * pal at on T, pivot oa 
ini a* tab to "W* aealo. Drw lino fro* "D" to "0" and pivot oa icalo "F" to point 
0 oa "J". Answor lo at intarieatioa of *FJ* w£t& *H*. 



Figure 10. - Homogram for calculation of &{.. (A. large working copy of this fig. rosy be obtained by 

using request card bound in back of report . ) 


UqolA t M y or* taro, 
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• 1 «2 #3 ft *5 *6 . 7 *8 *9 1,0 


Figure 11. - Construction of new axis j case of initial liquid temperature 
higher than 550° R for decane. Solid curve and new axis give the desired 
reduced curve. 
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Figure 12. - Uniform scales for construction of soales 10 Inches or less In length, 
copy of this fig. may be obtained by using request card bound In back of report.; 


(A large working 
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Draw line from T t 


scale. Pivot on intersection 


'tt' jc 

with. Q and draw line to point on 1-m . Intersection on r 
scale gives corresponding value of r*. 


t * Mass of drop at time 6 

l-m = ^ - .■7, — i ■ - - — 

Initial mass 

* _ Radius of drop at time 8 
r = Initial radius 

T l ( 0 ) = Initial temperature of liquid, 9 = 0, °R 
T L(0 ) = Temperature of liquid at time 9, °R 


1 -.2034 

r* 


Mil. 

r sod 


Mo) 

r 550 


Figure 13. - Nomogram for calculation of reduced radius from m , T w data. 



Percent mass transferred 



Time, sec 

fal Total urea sure . 1 atmosnhere absolute: a-t-r +.emrvovntnTv> . i 2ftn° r? »ir xn.1 nr.1 -Kr i c>n a a v^>, Q o 

»# — r ‘ / ~ ■-"*£ “■ } ■*</ - *■■■* ' Mj ; e.uv « W yUA 

second; initial radius, 0-108 inch; initial liquid temperature, 612° R, 


TMi 


1 A. _ HATOrko^ or«M a r\-P trn-»t-j ruin m<1 — .4-X am*. -LV i 1 J I -O a 

- wui^oj.j.uvwn wj- Tcujuww ucu.uuj.auxv.u£> wxou cApcx xmcju.ua, i. uaua xux - n— aecane- 
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X >XlX 




Experimental 

Reduced charts 

IBM calculation 

Solid sphere, steady state 

Steady state only 




Time , sec 


(h) Total pressure, 1 atmosphere absolute* air temperature, 1040° Bj air velocity, 97 inches per 

rtrifl ; ’T rv? +4 oT ffliHiia n ■< MftV> i 1 4 ^ ncO T) 


*j v • v^-v I AW4. <14.00. J-JU^lLLU UUlll^X aiyLLL C f kJ ( kJ JCV • 

Figure 14. - Continued. 
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totri pressure, i atmosphere absolute* air temperature, 1380° R; air velocity, 101.4- inches n 
second* initial radius, 0.0102 inch* initial liquid temperature, 650° R. 


Figure 14. - Concluded. 
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Figure 15. - Comparison of histories obtained, from (fig. 7) reduced-variable plot with IBM calculations. 
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(a) Air velocity, 900 e~®^; initial radius, 250 microns. 
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Figure 16. - Coraparison of reduced charts and IBM calculations for s-decane droplet with variable 
air velocity. Total pressure, 1 atmosphere absolute; air temperature, 1600° R; initial liquid 
temperature, 550° R. 
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Figure 17. - Flew diagram of IBM 650 program for vaporization calculations. 
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